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A theory describing the operation of a superconducting nanowire quantum interference device 
(NQUID) is presented. The device consists of a pair of thin-film superconducting leads connected 
by a pair of topologically parallel ultra-narrow superconducting wires. It exhibits intrinsic elec- 
trical resistance, due to thermally-activated dissipative fluctuations of the superconducting order 
parameter. Attention is given to the dependence of this resistance on the strength of an externally 
applied magnetic field aligned perpendicular to the leads, for lead dimensions such that there is 
essentially complete and uniform penetration of the leads by the magnetic field. This regime, in 
which at least one of the lead dimensions — length or width — lies between the superconducting co- 
herence and penetration lengths, is referred to as the mesoscopic regime. The magnetic field causes 
a pronounced oscillation of the device resistance, with a period not dominated by the Aharonov- 
Bohm effect through the area enclosed by the wires and the film edges but, rather, in terms of the 
geometry of the leads, in contrast to the well-known Little-Parks resistance of thin-walled supercon- 
ducting cylinders. A detailed theory, encompassing this phenomenology quantitatively, is developed 
through extensions, to the setting of parallel superconducting wires, of the Ivanchenko-Zil'berman- 
Ambegaokar-Halperin theory of intrinsic resistive fluctuations in a current-biased Josephson junc- 
tions and the Langer-Ambegaokar-McCumber-Halperin theory of intrinsic resistive fluctuations in 
superconducting wires. In particular, it is demonstrated that via the resistance of the NQUID, the 
wires act as a probe of spatial variations in the superconducting order parameter along the perimeter 
of each lead: in essence, a superconducting phase gradiometer. 

PACS numbers: 74.78.Na, 85.25.Dq, 85.65.+h 



I. INTRODUCTION 

The Little-Parks effect concerns the electrical resis- 
tance of a thin cylindrically-shaped superconducting film 
and, specifically, the dependence of this resistance on 
the magnetic flux threading the cylinder P, 0j- It is 
found that the resistance is a periodic function of the 
magnetic field, with period inversely proportional to the 
cross-sectional area of the cylinder. Similarly, in a DC 
SQUID, the critical value of the supercurrent is periodic 
in magnetic field, with period inversely proportional to 
the area enclosed by the SQUID ring 

In this Paper, we consider a mesoscopic analog of a 
DC SQUID. The analog consists of a device composed 
of a thin superconducting film patterned into two meso- 
scopic leads that are connected by a pair of (topologi- 
cally) parallel, short, weak, superconducting wires. Thus, 
we refer to the device as an NQUID (superconducting 
nanowire quantum interference device). The only restric- 
tion that we place on the wires of the device is that they 
be thin enough for the order parameter to be taken as 
constant over each cross-section of a wire, varying only 
along the wire length. In principle, this condition of one- 
dimensionality is satisfied if the wire is much thinner than 
the superconducting coherence length £. In practice, it is 
approximately satisfied provided the wire diameter d is 
smaller than 4.4 £ For thicker wires, vortices can ex- 
ist inside the wires, and such wires may not be assumed 
to be one dimensional. 

By the term mesoscopic we are characterizing phenom- 



ena that occur on length-scales larger than the supercon- 
ducting coherence length £ but smaller than the electro- 
magnetic penetration depth associated with magnetic 
fields applied perpendicular to the superconducting film. 
We shall call a lead mesoscopic if at least one of its two 
long dimensions is in the mesoscopic regime; the other di- 
mension may be either mesoscopic or macroscopic. Thus, 
a weak magnetic field applied perpendicular to a meso- 
scopic lead will penetrate the lead without appreciable 
attenuation and without driving the lead from the ho- 
mogeneous superconducting state to the Abrikosov vor- 
tex state. This is similar to the regime of operation of 
superconducting wire networks; see e.g., Ref. !■ The 
nanowires connecting the two leads are taken to be topo- 
logically parallel (i.e. parallel in the sense of electrical 
circuitry): these nanowires and edges of the leads define 
a closed geometrical contour, which will be referred to 
as the Aharonov-Bohm (AB) contour. In our approach, 
the nanowires are considered to be links sufficiently weak 
that any effects of the nanowires on the superconductiv- 
ity in the leads can be safely ignored. 

The theory presented here has been developed 
to explain experiments conducted on DNA-templated 
NQUIDs These experiments measure the electrical 
resistivity of a pair of superconducting nanowires sus- 
pended between long superconducting strips (see Fig.QJ. 
In them, a current source is used to pass DC current from 
a contact on the far end of the left lead to one on the far 
end of the right lead. The voltage between the contacts 
is measured (and the resistance is hence determined) as 
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FIG. 1: Schematic depiction of the superconducting phase 
gradiometer. A current J is passed through the bridges in 
the presence of a perpendicular magnetic field of strength B. 



a function of the magnetic field applied perpendicular to 
the plane of the strips. 

In the light of the foregoing remarks, the multiple- 
connectedness of the device suggests that one should an- 
ticipate oscillations with magnetic field, e.g., in the device 
resistance. Oscillations are indeed observed. But they 
are distinct from the resistance oscillations observed by 
Little and Parks and from the critical current oscillations 
observed in SQUID rings. What distinguishes the resis- 
tance oscillations reported in Ref. |5| from those found, 
e.g., by Little and Parks? First, the most notable aspect 
of these oscillations is the value of their period. In the 
Little-Parks type of experiment, the period is given by 
<&o/2ab, where $q(= hc/2e) is the superconducting flux 
quantum, 2a is the bridge separation, and b is the bridge 
length, i.e., the superconducting flux quantum divided 
by the area of the AB contour (see Fig. |2J). In a high- 
magnetic-field regime, such periodic behavior is indeed 
observed experimentally, with the length of the period 
somewhat shorter but of the same order of magnitude 
as in the AB effect ||. However, in a low-magnetic- 
field regime, the observed period is appreciably smaller 
(in fact by almost two orders of magnitude for our de- 
vice geometry). Second, because the resistance is caused 
by thermal phase fluctuations (i.e. phase slips) in very 
narrow wires, the oscillations are observable over a wide 
range of temperatures (~ IK). Third, the Little-Parks 
resistance is wholly ascribed to a rigid shift of the R(T) 
curve with magnetic field, as T c oscillates. In contrast, 
in our system we observe a periodic broadening of the 
transition (instead of the Little-Parks — type rigid shift) 
with magnetic field. Our theory explains quantitatively 
this broadening via the modulation of the barrier heights 
for phase slips of the superconducting order parameter 
in the nanowires. 

In the experiment, the sample is cooled in zero mag- 
netic field, and the field is then slowly increased while 
the resistance is measured. At a sample-dependent field 
(~ 5 mT) the behavior switches sharply from a low-field 
to a high-field regime. If the high-field regime is not 
reached before the magnetic field is swept back, the low- 
field resistance curve is reproduced. However, once the 
high-field regime has been reached, the sweeping back of 
the field reveals phase shifts and hysteresis in the R{B) 
curve. The experiments p] mainly address rectangular 
leads that have one mesoscopic and one macroscopic di- 



mension. Therefore, we shall concentrate on such strip 
geometries. We shall, however, also discuss how to ex- 
tend our approach to generic (mesoscopic) lead shapes. 
We note in passing that efficient numerical methods, such 
as the boundary element method (BEM) [|| , are available 
for solving the corresponding Laplace problems. 

This paper is arranged as follows. In Section [H] we 
construct a basic picture for the period of the magnetore- 
sistance oscillations of the two-wire device, which shows 
how the mesoscopic size of the leads accounts for the 
anomalously short magnetoresistance period in the low- 
field regime. In Section IlIII we concentrate on the prop- 
erties of mesoscopic leads with regard to their response 
to an applied magnetic field, and in Section II VI we ex- 
tend the LAMH model to take into account the inter-wire 
coupling through the leads. Analytical expressions are 
derived for the short- and long-wire limits, whilst a nu- 
merical procedure is described for the general case. The 
predictions of the model are compared with data from 
our experiment in Section and we give some conclud- 
ing remarks in Section lVll Certain technical components 
are relegated to the appendix, as is the analysis of exam- 
ple multiwire devices. 



II. ORIGIN OF MAGNETORESISTANCE 
OSCILLATIONS 

Before presenting a detailed development of the the- 
ory, we give an intuitive argument to account for the 
anomalously-short period of the magnetoresistance in the 
low-magnetic-field regime, mentioned above. 

A. Device geometry 

The geometry of the devices studied experimentally is 
shown in Fig. EI Five devices were successfully fabricated 
and measured. The dimensions of these devices are listed 
in Table HJ along with the short magnetoresistance oscil- 
lation period. The perpendicular penetration depth Aj_ 
for the films used to make the leads is roughly 70 /im, 
and coherence length £ is roughly 5 nm. 

B. Parametric control of the state of the wires by 
the leads 

The essential ingredients in our model are (i) leads, in 
which the applied magnetic field induces supercurrents 
and hence gradients in the phase of the order parame- 
ter, and (ii) the two wires, whose behavior is controlled 
parametrically by the leads through the boundary con- 
ditions imposed by the leads on the phase of the order 
parameters in the wires. For now, we assume that the 
wires have sufficiently small cross-sections that the cur- 
rents through them do not feed back on the order param- 
eter in the leads. (In Section llll Dl we shall discuss when 
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FIG. 2: Geometry of the two-wire device, showing the dimensions. The coordinate system used for the right lead (with the 
origin in the center of the lead) is also shown. The coordinates of the four corners of the right lead, as well as the coordinates 
of the points at which the two wires are connected to the right lead, are indicated. As shown, we always assume that the wires 
are attached near the center of the short edges of the leads. 
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FIG. 3: (a) Close-up of the two nanowires and the leads. The 
top (bottom) thick arrow represents the integration contour 
for determining the phase accumulation Oi^l^r (&2,l^r) in 
the first (second) wire. The dotted arrow in the left (right) 
lead indicates a possible choice of integration contour for de- 
termining the phase accumulation <52<-i, l (<52<-i, r)- These 
contours may be deformed without affecting the values of 
the various phase accumulations, as long as no vortices are 
crossed, (b) Sketch of the corresponding superconducting 
phase at different points along the AB contour when one vor- 
tex is located inside the contour. 



this assumption may be relaxed without altering the os- 
cillation period.) The dissipation results from thermally 
activated phase slips, which cause the superconducting 
order parameter to explore a discrete family of local min- 
ima of the free energy. (We assume that the barriers sep- 
arating these minima are sufficiently high to make them 
well-defined states.) These minima (and the saddle-point 
configurations connecting them) may be indexed by the 
net (i.e. forward minus reverse) number of phase slips 
that have occurred in each wire (n\ and n2, relative to 
some reference state). More usefully, they can be in- 
dexed by n s = min(ni, 712) (i.e. the net number of phase 
slips that have occurred in both wires) and n v = n\ — n-i 
(i.e. the number of vortices enclosed by the AB contour, 
which is formed by the wires and the edges of the leads). 
We note that two configurations with identical n v but 
distinct n s and n' s have identical order parameters, but 
differ in energy by 
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due to the work done by the current source supplying the 
current /, in which V is the inter-lead voltage, is the 
inter-lead phase difference as measured between the two 



points half-way between the wires, and the Josephson 
relation = 2eV/h has been invoked. In our model, we 
assume that the leads are completely rigid. Therefore the 
rate of phase change, and thus the voltage, is identical 
at all points inside one lead. For sufficiently short wires, 
n v has a unique value, as there are no stable states with 
any other number of vortices. 

Due to the screening currents in the left lead, induced 
by the applied magnetic field B (and independent of 
the wires), there is a field-dependent phase <$2<-i,l(-B) = 

J. dr ■ Vip(B) (computed below) accumulated in pass- 
ing from the point at which wire 1 (the top wire) con- 
tacts the left (L) lead to the point at which wire 2 (the 
bottom wire) contacts the left lead (see Fig. Simi- 
larly, the field creates a phase accumulation 52<—i,r(B) 
between the contact points in the right (R) lead. As the 
leads are taken to be geometrically identical, the phase 
accumulations in them differ in sign only: 82^- i,l{B) = 
—^2<-i,r(B). We introduce 6(B) = 62^1. l(B). In deter- 
mining the local free-energy minima of the wires, we solve 
the Ginzburg-Landau equation for the wires for each vor- 
tex number n v , imposing the single- valuedness condition 
on the order parameter, 
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h.L^R + 25(B) = 2nn v 



(2) 



This condition will be referred to as the phase constraint. 
Here, 6\.l^r = f^dr- Vip(B) is the phase accumulated 
along wire 1 in passing from the right to the left lead; 
02.l^r. is similarly defined for wire 2. 

Absent any constraints, the lowest energy configura- 
tion of the nanowires is the one with no current through 
the wires. Here, we adopt the gauge in which A = Bye x 
for the electromagnetic vector potential, where the coor- 
dinates are as shown in Fig. [21 The Ginzburg-Landau 
expression for the current density in a superconductor is 
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For our choice of gauge, the vector potential is always 
parallel to the nanowires, and therefore the lowest en- 
ergy state of the nanowires corresponds to a phase ac- 
cumulation given by the flux through the AB contour, 
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&i,l^b = —0 2 ,l^r — 2irBab/$ - As we shall show 
shortly, for our device geometry (i.e. when the wires are 
sufficiently short, i.e., b <C 0; this phase accumulation 
may be safely ignored, compared to the phase accumu- 
lation 8(B) associated with screening currents induced 
in the leads. As the nanowires are assumed to be weak 
compared to the leads, to satisfy the phase constraint J3J, 
the phase accumulations in the nanowires will typically 
deviate from their optimal value, generating a circulat- 
ing current around the AB contour. As a consequence 
of LAMH theory, this circulating current results in a de- 
crease of the barrier heights for phase slips, and hence an 
increase in resistance. The period of the observed oscilla- 
tions is derived from the fact that whenever the magnetic 
field satisfies the relation 
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[where m is an integer and the factor of 2 accompanying 
5(B) reflects the presence of two leads], there is no cir- 
culating current in the lowest in energy state, resulting 
in minimal resistance. Furthermore, the family of free 
energy-minima of the two-wire system (all of which, in 
thermal equilibrium, are statistically populated accord- 
ing to their energies) is identical to the B = case. The 
mapping between configurations at zero and nonzero B 
fields is established by a shift of the index n v — > n v — m. 
Therefore, as the sets of physical states of the wires are 
identical whenever the periodicity condition |0J is sat- 
isfied, at such values of B the resistance returns to its 
B = value. 



C. Simple estimate of the oscillation period 

In this subsection, we will give a "back of the envelope" 
estimate for the phase gain 8(B) in a lead by consider- 
ing the current and phase profiles in one such lead. Ac- 
cording to the Ginzburg-Landau theory, in a mesoscopic 
superconductor, subjected to a weak magnetic field, the 
current density is given by Eq. J3J|. Now consider an 
isolated strip-shaped lead used in the device. Far from 
either of the short edges of this lead, A = Bye x is a 
London gauge 0, i.e., along all surfaces of the supercon- 
ductor A is parallel to them; A — > in the center of the 
superconductor; and V • A = 0. In this special case, the 
London relation [20| states that the supercurrent den- 
sity is proportional to the vector potential in the London 
gauge. Using this relation, we find that the supercurrent 
density is J oc — (2e/h) A = ~(2e/h) Bye x , i.e., there is 
a supercurrent density of magnitude oc (2e/H) Bl flowing 
to the left at the top (long) edge of the strip and to the 
right at the bottom (long) edge. At the two short ends 
of the strip, the two supercurrents must be connected, so 
there is a supercurrent density of magnitude ~ (2e/h)Bl 
flowing down the left (short) edge of the strip and up the 
right (short) edge (see Fig. @}. Near the short ends of 
the strips, our choice of gauge no longer satisfies the cri- 
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FIG. 4: Current profile in a long superconducting strip, cal- 
culated for a finite-length strip by summing the series for 
^4> ~ ^r~A [from Eq. I|15|l ] numerically. Note that there is 
no vortex in the center of the lead. 



tcria for being a London gauge, and therefore may be 
nonzero. As, in our choice of gauge, A points in the e x 
direction, the supercurrent on the ends of the strip along 
e y must come from the V y 4> term. Near the center of the 
short edge V y <f) = —2ircil/$oB. The phase difference be- 
tween the points (— L, —a) and (— L, a) is therefore given 
by 



V y (t>dy = — ^-B2al, 

-a *0 



(5) 



where we have substituted 2-7r/$o for 2e/h and C\(a/l) 
is a function of order unity, which accounts for how the 
current flows around the corners. As we shall show, c\ 
depends only weakly on a//, and is constant in the limit 
a < I. 

Finally, we obtain the magnetoresistance period by 
substituting Eq. JSJ into Eq. (QJ: 
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Thus, we see that for certain geometries the period is 
largely determined not by the flux threading through 
the geometric area 2ab but by the response of the leads 
and the corresponding effective area Aal, provided the 
nanowires are sufficiently short (i.e. & < !), justifying 
our assumption of ignoring the phase gradient induced 
in the nanowires by the magnetic field. 

In fact, we can also make a prediction for the peri- 
odicity of the magnetoresistance at high magnetic fields, 
i.e., when vortices have penetrated the leads (see Sec- 
tion ^^^J. To do this, we should replace I in Eq. © 
by the characteristic inter- vortex spacing r. Note that 
if r is comparable to 6, we can no longer ignore the flux 
through the AB contour. Furthermore, if r -C b then 
the flux through the AB contour determines periodicity 
and one recovers the usual Aharonov-Bohm type of phe- 
nomenology. 
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III. MESOSCALE SUPERCONDUCTING 
LEADS 

In this section and the following one we shall develop a 
detailed model of the leads and nanowires that constitute 
the mesoscopic device. 



slow oscillations is found to occur at 3.1 mT for that 
sample 5j. It is possible to determine the critical mag- 
netic fields H s and H c i by the direct imaging of vor- 
tices. Although we do not know of such a direct measure- 
ment of H Sl H c \ was determined by field cooling niobium 
strips, and found to agree in magnitude to Likharev's es- 
timate [lCj . 



A. Vortex- free and vorticial regimes 

Two distinct regimes of magnetic field are expected, 
depending on whether or not there are trapped (i.e. lo- 
cally stable) vortices inside the leads. As described by 
Likharev |8(, a vortex inside a superconducting strip- 
shaped lead is subject to two forces. First, due to the the 
currents induced by the magnetic field there is a Magnus 
force pushing it towards the middle of the strip. Second, 
there is a force due to image vortices (which are required 
to enforce the boundary condition that no current flows 
out of the strip and into the vacuum) pulling the vortex 
towards the edge. When the two forces balance at the 
edge of the strip, there is no energy barrier preventing 
vortex penetration and vortices enter. Likharev has esti- 
mated of the corresponding critical magnetic field to be 



nd £a(l) 



(7) 



where d(= 21) is the width of the strip and a(l) ~ 1 
for strips that are much narrower than the penetration 
depth (i.e. for d<A). 

Likharev has also shown that, once inside a strip, vor- 
tices remain stable inside it down to a much lower mag- 
netic field H c i, given by 



H <* = — 73 m 77 
nd d \4£ 



(8) 



At fields above H c i the potential energy of a vortex in- 
side the strip is lower than for one outside (i.e. for a 
virtual vortex 0). Therefore, for magnetic fields in the 
range H c \ < H < H s vortices would remain trapped in- 
side the strip, but only if at some previous time the field 
were larger than H s . This indicates that hysteresis with 
respect to magnetic field variations should be observed, 
once H exceeds H s and vortices become trapped in the 
leads. 

In real samples, in addition to the effects analyzed by 
Likharev, there are also likely to be locations (e.g. struc- 
tural defects) that can pin vortices, even for fields smaller 
than H c i , so the reproducibility of the resistance vs. field 
curve is not generally expected once H s has been sur- 
passed. 

As magnetic field at which vortices first enter the leads 
is sensitive to the properties of their edges, we expect only 
rough agreement with Likharev's theory. For sample 219- 
4, using Likharev's formula, we estimate H s = 11 mT 
(with £ = 5nm). The change in regime from fast to 



B. Phase variation along the edge of the lead 

In the previous section it was shown that the period- 
icity of the magnetoresistance is due to the phase accu- 
mulations associated with the currents along the edges of 
the leads between the nanowires. Thus, we should make 
a precise calculation of the dependence of these currents 
on the magnetic field, and this we now do. 



1. Ginzburg-Landau theory 

To compute 5(B), we start with the Ginzburg-Landau 
equation for a thin film as our description of the meso- 
scopic superconducting leads: 

2 



atp + /3| VI W - 
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ip = 0. (9) 



Here, if) is the Ginzburg-Landau order parameter, e* (= 
2e) is the charge of a Cooper pair and m* is its mass, 
and a and /? may be expressed in terms of the coherence 
length £ and critical field H c via a = — ti 2 /2m*£ 2 and 
f3 = 4Tra 2 /Hl 

The assumptions that the magnetic field is sufficiently 
weak and that the lead is a narrow strip (compared with 
the magnetic penetration depth) allow us to take the 
amplitude of the order parameter in the leads to have 
the value appropriate to an infinite thin film in the ab- 
sence of the field. By expressing the order parameter in 
terms of the (constant) amplitude t/>o and the (position- 
dependent) phase 4>(r), i.e., 

V(r)=Voe^W, (10) 
the Ginzburg-Landau formula for the current density, 
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-r^A(r), (11) 



becomes 



J = ±-4, 2 (hV^r)~^-A(r)), (12) 

and [after dividing by e l ^ r ^] the real and imaginary parts 
of the Ginzburg-Landau equation become 
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(13a) 
(13b) 
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As long as any spatial inhomogeneity in the gauge- 
covariant derivative of the phase is weak on the length- 
scale of the coherence length [ i.e. flV^(r) — j^A(r) <C 
1 ] , the third term in Eq. (jl3a[l is much smaller than the 
first two and may be ignored, fixing the amplitude of the 
order parameter at its field-free infinite thin film value, 
viz., -0o = y/~ a//3- To compute 4>(r) we need to solve 
the imaginary part of the Ginzburg-Landau equation. 




2. Formulation as a Laplace problem 

We continue to work in the approximation that the 
amplitude of the order parameter is fixed at tpQ. Start- 
ing from Eq. I|13b(l . we see that for our choice of gauge, 
A = Bye x , the phase of the order parameter satisfies 
the Laplace equation, V 2 = 0. We also enforce the 
boundary condition that no current flows out of the su- 
perconductor on boundary surface E, whose normal is n: 
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FIG. 5: Phase profile in the leads in the vicinity of the trench, 
generated by numerically summing the series for <j> for a finite- 
length strip. Arrows indicate phases connected by nanowires. 



Eq. 1|15|) taking the values 

h (-1)" 



B k = -A k = 



k^l cosh(fc„L 

and hence to the solution 

(-l)"2/i 
k^l cosh(fc n L) 



.») = £ 



(71 = 0,1,...), (16) 



sm(k n y) smh(k n x). (17) 



3. Solving the Laplace problem for the strip geometry 



Figure |S] shows the phase profiles in the leads, in the 
region close to the trench that separates the leads. 



To solidify the intuition gained via the physical argu- 
ments given in Section [n] we now determine the phase 
profile for an isolated superconducting strip in a magnetic 
field. This will allow us to determine the constant c\ in 
Eq. © , and hence obtain a precise formula for the mag- 
netoresistance period. To this end, we solve Laplace's 
equation for <f> subject to the boundary conditions (|14fl . 
We specialize to the case of a rectangular strip plf . 

In terms of the coordinates defined in Fig-EI we expand 
<p(x, y) as the superposition 

4>{x, y) = 9 L/R + Y, ( A k e~ kx + B k e kx ) sin(fcy), (15) 



which automatically satisfies Laplace's equation, al- 
though the boundary conditions remain to be satisfied. 
©l(r) is the phase at the the point in the left (right) 
lead located half-way between the wires. In other words 
Ql = <P(—L — b, 0) and Or = <p{—L, 0) in the coordinate 
system indicated in Fig. 6l/r & r e not determined by 
the Laplace equation and boundary conditions, but will 
be determined later by the state of the nanowires. 

We continue working in the gauge A = By e x . The 
boundary conditions across the edges at y = ±1 (i.e. the 
long edges) are d y 4>(x,y = ±1) = 0. These conditions 
are satisfied by enforcing k n — ir(n + where n = 

0, 1,2,.. .. The boundary conditions across the edges at 
x = ±L (i.e. the short edges) are d x <p(x = ±L,y) = hy 
(where h = 2ttB/^q). This leads to the coefficients in 



C. Period of magnetoresistance for leads having a 
rectangular strip geometry 

Using the result for the phase that we have just estab- 
lished, we see that the phase profile on the short edge of 
the strip at x = — L is given by 

T . 2^ 2 ^ (-1)" . 7r(n+i)y 
${-L,y) = an , . ( 18 

where we have taken the limit L — > oo. We would like to 
evaluate this sum at the points (x, y) = (— L, ±a). This 
can be done numerically. For nanowires that are close 
to each other (i.e. for a <^ I), an approximate value can 
be found analytically by expanding in a power series in 
a around y — 0: 



d 

-L, a) = cj)(-L, 0) + a —tft-L, y) 



y=0 



,2 Q2 



(19) 



+ 0(a 3 ). 



The first and third terms are evidently zero, as <f> is an 
odd function of y. The second term can be evaluated 
by changing the order of summation and differentiation. 
(Higher-order terms are harder to evaluate, as the chang- 
ing of the order of summation and differentiation does not 
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FIG. 6: Phase profile on the x = —L (i.e. short) edge of the 
strip. Numerical summation [Eq. I|18|l with 100 terms] for 
L — 2 and L = oo (I = 1) as well as the linear form from 
Eq. 12011 . Note that the linear fit is good near the origin 
(e.g. for a < 0.25Z), and the curves for L = 2 and L = oo 
coincide. 



work for them.) Thus, to leading order in a we have 

8G 



(-L,o) 



-hla. 



(20) 



where G = V°° n ^ ~ 0.916 is the Catalan num- 

ber (see Ref. ^2). This linear approximation is plotted, 
together with the actual phase profile obtained by the 
numerical evaluation of Eq. (|18fl . in Fig. Hence, the 
value of ci in Eq. JSJ becomes C\ — 8G/n 2 w 0.74, and 
Eq. jnjl becomes 



2^' 



7T 2 $0 

8G 4aZ ' 



(21) 



To obtain this result we used the relation S(B)/2 = 
<f>(-L,a). 



d x cj) = /($ /-ff c 2 s£)(<5(y + a) - S(y - a)), where / is the 
current circulating in the loop, and H c , s, £ are the film 
critical field, thickness, and coherence length. By using 
the same expansion as before, Eq. (|15|) . we obtain the 
coefficients of the Fourier-series in the long strip limit: 



A k = I(<l>o/Hls(;) 



2 sin(fca) 
kl exp(fcZ) 



(22) 



Having the coefficients of the Fourier series, we can find 
the phase difference in the right lead between the two 
points at which the nanowires connect to the right lead, 
induced in this lead by the current circulating in the loop: 



S cc = 2/($o/# c V) 



k=l/w 

E 

n=0 



sin 2 (fca) 
kl 



ln(2//7rw). (23) 



Here, we have introduced a large wave- vector k cut-off 
at the inverse of the width w of the wire. On the other 
hand, the current flowing through the wire is 



A9 



$ b 



(24) 



where £, H c , and s are the wire coherence length, critical 
field, and height (recall that b is the wire length). To 
support a circulating current that corresponds to a phase 
accumulation of AO along one of the wires, the phase 
difference between the two nanowires in the lead must be 
on the order of 



Sec. = A8 



b (-HM)fi 



film 



7TW 



(25) 



For our experiments , we estimate that the ratio of S cc 
to Ad is always less than 20%, validating the assumption 
of weak coupling. 



D. Bridge- lead coupling 

In order to simplify our analysis we have assumed that 
the nanowires do not exert any influence on the order 
parameter in the leads. We examine the justification for 
this assumption in the setting of the experiment that we 
are attempting to describe (3- 

The assumption will be valid if the bending of the 
phase of the order parameter, in order to accommodate 
any circulating current around the AB contour, occurs 
largely in the nanowires. As the phase of the order pa- 
rameter in the leads satisfies the Laplace equation, which 
is linear, we can superpose the circulating-current solu- 
tion with the previously-obtained magnetic-field-induced 
solution. The boundary conditions on the right lead for 
the circulating-current solution are n ■ V(/> = every- 
where, except at the two points where the nanowires are 
attached to the lead [i.e. at (x,y) — (L,±a)]. Treat- 
ing the nanowires as point current sources, the bound- 
ary condition on the short edge of the right lead is 



E. Strong nanowires 

We remark that the assumption of weak nanowires is 
not obligatory for the computation the magnetoresistance 
period. Dropping this assumption would leave the period 
of the magnetoresistance oscillations unchanged. 

To see this, consider (j>n, i.e., the phase profile in the 
leads that corresponds to the lowest energy solution of 
the Ginzburg-Landau equation at field corresponding to 
the first resistance minimum [i.e. at B being the first 
non-zero solution of Eq. |@}]. For this case, and for 
short wires, the phase gain along the wires is negligi- 
ble, whereas the phase gain in the leads is 2tt 7 even for 
wires with large critical current. Excited states, with 
vortices threading the AB contour, can be constructed 
by the linear superposition of 4>n with </>on„, where </>on„ 
is the phase profile with n v vortices at no applied mag- 
netic field. 

This construction requires that the nanowires are nar- 
row, but works independently of whether nanowires are 
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strong or weak, in the limit that H <C H c . The energy of 
the lowest energy state always reaches its minimum when 
the applied magnetic field is such that there is no phase 
gain (i.e. no current) in the nanowires. By the above 
construction, it is clear that the resistance of the device 
at this field is the same as at zero field, and therefore the 
minimum possible. 

Therefore, our calculation of the period is valid, inde- 
pendent of whether the nanowires are weak or strong. 
However, the assumption of weak nanowires is neces- 
sary for the computation of magnetoresistance ampli- 
tude, which we present in the following section. 

IV. PARALLEL SUPERCONDUCTING 
NANOWIRES AND INTRINSIC RESISTANCE 

In this section we consider the intrinsic resistance of 
the device. We assume that this resistance is due to 
thermally activated phase slips (TAPS) of the order pa- 
rameter, and that these occur within the nanowires. 
Equivalently, these processes may be thought of as 
thermally activated vortex flow across the nanowires. 
Specifically, we shall derive analytical results for the 
asymptotic cases of nanowires that are either short 
or long, compared to coherence length, i.e. Josephson 
junctions |l2l Il3| or Langer-Ambegaokar-McCumber- 
Halperin (LAMH) wires [3 [13 ; see also Ref. O- We 
have not been able to find a closed-form expression for 
the intrinsic resistance in the intermediate-length regime, 
so we shall consider that case numerically. 

There are two (limiting) kinds of experiments that may 
be performed: fixed total current and fixed voltage. In 
the first kind, a specified current is driven through the 
device and the time-averaged voltage is measured. Here, 
this voltage is proportional to the net number of phase 
slips (in the forward direction) per unit time, which de- 
pends on the height of the free-energy barriers for phase 
slips. Why do we expect minima in the resistance at 
magnetic fields corresponding to 26 = 2rmr and max- 
ima at 25 = (2m + l)ir for m integral, at least at van- 
ishingly small total current through both wires? For 
25 = 2to7t the nanowires are unfrustrated, in the sense 
that there is no current through either wire in the lowest 
local minimum of the free energy. On the other hand, 
for 26 = (2m + 1)tt the nanowires are maximally frus- 
trated: there is a nonzero circulating current around the 
AB contour. Quite generally, the heights of the free- 
energy barriers protecting locally stable states decrease 
with increasing current through a wire, and thus the frus- 



trated situation is more susceptible to dissipative fluctua- 
tions, and hence shows higher resistance. Note, however, 
that due to the inter-bridge coupling caused by the phase 
constraint, the resistance of the full device is more subtle 
than the mere addition of the resistances of two inde- 
pendent, parallel nanowires, both carrying the requisite 
circulating current. 

In the second kind of experiment, a fixed voltage is 
applied across the device and the total current is mea- 
sured. In this situation, the inter-lead voltage is fixed, 
and therefore the phase drop along each wire is a fixed 
function of time. Hence, there is no inter-bridge coupling 
in the fixed voltage regime. Therefore, the resistance of 
the device would not exhibit magnetic field dependence. 
If the voltage is fixed far away from the wires, but not 
in the immediate vicinity of the wires, so that the phase 
drop along each wire is not rigidly fixed, then some of 
the magnetic field dependence of the resistance would be 
restored. In our experiments on two-wire devices, we be- 
lieve that the situation lies closer to the fixed current 
limit than to the fixed voltage limit, and therefore we 
shall restrict our attention to the former limit. 

In the fixed-current regime, the relevant independent 
thermodynamic variable for the device is the total current 
through the pair of wires, i.e., / = I\+l2- Therefore, the 
appropriate free energy to use, in obtaining the barrier 
heights for phase slips, is the Gibbs free energy G(I), as 
discussed by McCumber ^tJi rather than the Helmholtz 
free energy F(0) j^. In the Helmholtz free energy the 
independent variable can be taken to be = ©l — 0r, 
i.e., the phase difference across the center of the "trench," 
defined modulo 2tt. G(I) is obtained from F{&) via the 
appropriate Legendre transformation: 

G(I)=F(Q)- A/e, (26) 

where the second term represents the work done on the 
system by the external current source. F(<3) is the sum of 
the Helmholtz free energies for the individual nanowires: 



F(Q) = Fi(0i) + F 2 (0 2 ), (27) 

where -Fi(2)($i(2)) is the Ginzburg-Landau free energy for 
first (second) wire and a simplified notation has been 
used 6*i = 9i l^r and 62 = 6*2, l^b- #1 and 62 are re- 
lated to each other and to through the phase con- 
straint Eq. J2J. 



I 

A. Short nanowires: Josephson junction limit 

If the nanowires are sufficiently short, they may be treated as Josephson junctions. Unlike the case of long nanowires, 
described in the following subsection, in this Josephson regime there is no metastability, i.e., the free energy of each 
junction is a single-valued function of the phase difference, modulo 2-7r, across it. The phase constraint then implies 
that there is a rigid difference between the phases across the two junctions. As a consequence, n v can be set to zero. 
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T (K) B(uT) 

FIG. 7: Sample 219-4: Experimental data (solid lines) and theoretical fits using the Josephson junction model (dashed lines), 
(a): Resistance vs. temperature curves. (1) Zero magnetic field and low total current. (2) Magnetic field set to maximize the 
resistance and low total current. (3) Zero magnetic field and 70 nA total current, (b): Resistance as a function of magnetic field 
at various temperatures from 1.2 to 2.0 K in 0.1 K increments. The fitting parameters used were J c i = 639 nA, J C 2 = 330 nA, 
Tci = 2.98 K, and T C 2 = 2.00 K, with corresponding coherence lengths £i(0) = 23 nm and £2(0) = 30 nm. Only one set of fitting 
parameters [derived from curves (1) and (2)] was used to produce all the theoretical curves. 



The Gibbs free energy in such a configuration is then 

h 



G(I) = ~Y e { IcX cos(0l) + /c2 cos(02) + /9 ) ' (28) 

where I c i and I c2 are the critical currents for the junctions. In thermodynamic equilibrium, the Gibbs free energy 
must be minimized, so the dependent variable 8 must be chosen such that dG(I)/d<d = 0. 
Using 61 = Q + 5 and 82 = — 6, G(I) may be rewritten in the form 



G(I) = + Ic2) 2 co S 2 S + (I cl - I c2 ) 2 sin 2 S ■ cos(tf) (29) 



/cl-Jc2 
/c2+/cl 



tan 6 



and ■& = 6 



where we have shifted the free energy by an additive constant (^)Jtan 1 

tan -1 ( 7^+7^ ) tan J . In this model, the option for having I c \ 7^ J C 2 is kept open. Equation (|29|l shows that, 

up to an additive constant, the free energy of the two-junction device is identical to that of an effective single-junction 
device with an effective I Cl which is given by 



I c2 ) 2 cos 2 6 + (id - I c2 ) 2 sin" 6. (30) 

Thus, we may determine the resistance of the two-junction device by applying the well-known results for a single 
junction, established by Ivanchenko and Zil'berman |l2j and by Ambegaokar and Halperin fl3| : 

2(1 -x 2 ) 1 / 2 / / 1 \ 

R = R n — exp ^-7(V 1 - x 2 + x skT 1 x)J sinh(7T7x/2), (31a) 

x = I/I c , 7 = hI c /ek B T, (31b) 

where R n is the normal-state resistance of the two-junction device. This formula for R holds when the free-energy 
barrier is much larger than ksT, so that the barriers for phase slips are high. References [l2l IT^ | provide details 
on how to calculate the resistance in the general case of an over damped junction, which includes that of shallow 
barriers. Figurc|7|shows the fits to the resistance, computed using Eqs. H30II31|) . as a function of temperature, magnetic 
field, and total current for sample 219-4. Observe that both the field- and the temperature-dependence are in good 
agreement with experimental data. In Section IV B 21 we make more precise contact between theory and experiment, 
and explain how the data have been fitted. We also note that, as it should, our Josephson junction model exactly 
coincides with our extension of the LAMH model in the limit of very short wires and for temperatures for which the 
barrier-crossing approximation is valid. 

B. Longer nanowires: LAMH regime Il5| . which we shall use to make a quantitative estimate 

In this section we describe an extension of the LAMH 
model of resistive fluctuations in a single narrow wire [lH 
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of the voltage across the two- wire device at a fixed total 
current. In this regime the nanowires are sufficiently long 
that they behave as LAMH wires. We shall only dwell on 
two-wire systems, but we note in passing that the model 
can straightforwardly be extended to more complicated 
sets of lead interconnections, including periodic, grating- 
like arrays (see Appendix 0. 

As the sample is not simply connected, i.e., there is 
a hole inside the AB contour, it is possible that there 
are multiple metastable states that can support the total 
current. These states differ by the number of times the 
phase winds along paths around the AB contour. The 
winding number n v changes whenever a vortex (or an 
anti- vortex) passes across one of the wires. 

In the present theory, we include two kinds of processes 
that lead to the generation of a voltage difference between 
the the leads; see Fig. [5] In the first kind of process 
(Fig. |Ht), two phase slips occur simultaneously: a vortex 
passes across the top wire and, concurrently, an anti- 
vortex passes across the bottom wire (in the opposite di- 
rection), so that the winding number remains unchanged. 
In the second kind of process (Fig.ISt*), the phase slips oc- 
cur sequentially: a vortex (or anti-vortex) enters the AB 
contour by passing across the top (or bottom) wire, stays 
inside the contour for some time- interval, and then leaves 
the AB contour through the bottom (or top) wire 

Our goal is to extend LAMH theory to take into ac- 
count the influence of the wires on each other. In Ap- 
pendix El we review some necessary ingredients associ- 
ated with the LAMH theory of a single wire. As the wires 
used in the experiments are relatively short (i.e. 10 to 20 
zero-temperature coherence lengths in length), we also 
take care to correctly treat the wires as being of finite 
length. 

Recall that we are considering experiments performed 
at a fixed total current, and accordingly, in all configura- 
tions of the order parameter this current must be shared 
between the top and bottom wires. We shall refer to this 
sharing, 

I = h+h, (32) 

as the total current constraint. Let us begin by consider- 
ing a phase-slip event in a device with an isolated wire. 
While the order parameter in that wire pinches down, 
the end-to-end phase accumulation must adjust to main- 
tain the prescribed value of the current through the wire. 
Now consider the two-wire device, and consider a phase 
slip event in one of the wires. As in the single-wire case, 
the phase accumulation will adjust, but in so doing it 
will alter the current flowing through other wire. Thus, 
in the saddle-point configuration of the two-wire system 
the current splitting will differ from that in the locally 
stable initial (and final) state. 

Taking into account the two kinds of phase-slip pro- 
cesses, and imposing the appropriate constraints (i.e. the 
total current constraint and the phase constraint), we 
construct the possible metastable and saddle-point con- 
figurations of the order parameter in the two- wire system. 



Finally, we compute the relevant rates of thermally ac- 
tivated transitions between these metastable states, con- 
struct a Markov chain 0], and determine the steady- 
state populations of these states. Thus, we are able to 
evaluate the time-average of the voltage generated be- 
tween the leads at fixed current due to these various 
dissipative fluctuations. We mention that we have not 
allowed for wires of distinct length or constitution (so 
that the Ginzburg-Landau parameters describing them 
are taken to be identical) . This is done solely to simplify 
the analysis; extensions to more general cases would be 
straightforward but tedious. 



1. Parallel pair of nanowires 

The total Gibbs free energy for the two- wire system is 
given by 

G(I) = Fi (6>i ) + F 2 (9 2 ) - 4£6 • [J x + J 2 ). (33) 

Here, we have followed MH by rewriting the current- 
phase term in terms of dimensionless currents in wires 
i = 1,2, i.e., Ji defined via Ii — 8ircJi£ /<I>o. More- 
over, E = — f^- is the condensate energy density per unit 
length of wire, and Fi{9i) is the Helmholtz free energy 
for a single wire along which there is a total phase ac- 
cumulation of 9i. The precise form of Fi(9i) depends on 
whether the wire is in a metastable or saddle-point state. 

We are concerned with making stationary the total 
Gibbs free energy at specified total current /, subject 
to the phase constraint, Eq. J2J). This can be accom- 
plished by making stationary the Helmholtz free energy 
on each wire, subject to both the total current constraint 
and the phase constraint, but allowing B\ and 9 2 to vary 
so as to satisfy these constraints — in effect, adopting the 
total current / as the independent variable. The station- 
ary points of the Helmholtz free energy for a single wire 
are reviewed in Appendix [U] as implicit functions of (9; , 
i.e., the end-to-end phase accumulation along the wire. 
The explicit variable used there is Ji, which is related to 
6i via Eq. JC2T)) . 



2. Analytical treatment in the limit of long nanowires 

In the long- wire limit, we can compute the resistance 
analytically by making use of the single-wire free energy 
and end-to-end phase accumulation derived by L anger 
and Ambegaokar 01 (and extended by McCumber [l7| 
for the case of the constant-current ensemble). Through- 
out the present subsection we shall be making an expan- 
sion in powers of 1/6, where 6 is the length of the wire 
measured in units of the coherence length, keeping terms 
only to first order in 1/6. Thus, one arrives at formulae for 
the end-to-end phase accumulations and Helmholtz free 
energies for single-wire metastable (m) and saddle-point 
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FIG. 8: Thermally activated phase slip processes under consideration, (a) Parallel phase slips, (b) Sequential phase slips. 



(sp) states QJ]: 
m («) = Kb, 



_i / 1 — 3k' 
/ sp {K) = kO + 2 tan [ 2 

F ro (/s) = -£ (6(1 - K 2 ) 2 ) , 
Fsp(«) = -£ (b(l-K 2 ) 2 -^ 



2\ 1/2 



VI - 3k 2 



(34a) 

(34b) 
(34c) 

(34d) 



where K is defined via Jj = /tj(l — k 2 ). In the small- 
current limit, one can make the further simplification 
that J, w k,;; henceforth we shall keep terms only up 
to first order in k. To this order, the phase difference 
along a wire in a saddle-point state becomes 



Kb + it - 2V2k. 



(35) 



Next, we make use of these single- wire LAMH results 
to find the metastable and saddle-point states of the two- 
wire system, and use them to compute the correspond- 
ing barrier heights and, hence, transition rates. At low 
temperatures, it is reasonable to expect that only the 
lowest few metastable states will be appreciably occu- 
pied. These metastable states, as well as the saddle-point 
states between them, correspond to pairs, k± and K2, one 
for each wire, that satisfy the total current constraint as 
well as the phase constraint: 



K\ + k 2 = J, 
h(*i) ~ 2 {k 2 ) = 2nn v + 25, 



(36) 
(37) 



where we need to substitute the appropriate O m / sp {Ki) 
from Eqs. l|51all3l)fl for 6^). 

In the absence of a magnetic field (i.e. 5 = 0), the 
lowest energy state is the one with no circulating current, 
and the current split evenly between the two wires. This 
corresponds to the solution of Eqs. I|3t)l I37[) with n = 0, 
together with the substitution l|34a|) for 0i(/tj) for both 
wires (i.e. 9\ = K\b and 02 = K 2 b). Thus we arrive at the 
solution: 



Ki = J/2, 
k 2 = J/2, 



0! = bJ/2, (38a) 
2 = bJ/2. (38b) 



If we ignore the lowest (excited) metastable states then 
only a parallel phase-slip process is allowed. The saddle 
point for a parallel phase slip corresponds to a solution 



of Eqs. with n = and the substitution (gSJl for 

9i(Ki) for both wires: 

Ki = J/2, 0! = bJ/2 + n- 2V2J/2, (39a) 
k 2 = J/2, 2 = bJ/2 + 7T - 2V2J/2. (39b) 

The change in the phase difference across the center of 
the trench, AO = [O sp — O m ], is 7r — 2y2K for a forward 
phase slip, and — ir — 2\/2k for a reverse phase slip. The 
Gibbs free-energy barrier for the two kinds of phase slips, 
computed by subtracting the Gibbs free energy for the 
ground state from that of the saddle-point state, is 



AG = £ 



16V2 



±4Jtt 



(40) 



The former free-energy is obtained by substituting 
Eq. H34c|l into Eq. I|33[l for both wires; the latter one is 
obtained by substituting Eq. (|34d|) into Eq. (|3^|) for both 
wires. We note that the Gibbs free-energy barrier heights 
for parallel phase slips (in both the forward and reverse 
directions) are just double those of the LAMH result for 
a single wire. From the barrier heights, we can work out 
the generated voltage by appealing to the Josephson re- 
lation, V — (H/2e)Q, and to the fact that each phase slip 
corresponds to the addition (or subtraction) of 27r to the 
phase. Hence, we arrive at the current-voltage relation 
associated with parallel phase slips at S = 0: 



, par 



= ^ e -^^ s inh(/// ), (41) 
e 



where the prefactor f2 may be computed using time- 
dependent Ginzburg-Landau theory or extracted from 
experiment, and Io = Ae/ (3h. 

If we take into account the two lowest excited states, 
which we ignored earlier, then voltage can also be gener- 
ated via sequential phase slips (in addition to the paral- 
lel ones, treated above). To tackle this case, we con- 
struct a diagram in which the vertices represent the 
metastable and saddle-point solutions of Eqs. I|3til I37fl . 
and the edges represent the corresponding free energy 
barriers; see Fig. Pairs of metastable-state vertices 
are connected via two saddle-point-state vertices, corre- 
sponding to a phase slip on either the top or the bottom 
wire. To go from one metastable state to another, the sys- 
tem must follow the edge out of the starting metastable 
state leading to the desired saddle-point state. We as- 
sume that, once the saddle-point state is reached, the 
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FIG. 9: Diagram representing the ground state (central filled circle), the two lowest-energy metastable states (left and right filled 
circles), and the saddle-point states connecting them (open circles) for the case of magnetoresistance minimum (i.e. 5 — 0, tt, . . .). 
The saddle-point states at the top of the graph correspond to phase slips in the top wire (i.e. wire 1); those at the bottom 
correspond to the bottom wire (i.e. wire 2). (Saddle-point states for parallel phase slips are not shown.) For each state Ki 
and K2 are listed. For each barrier (represented by an edge and labeled by a through h), the table at the right lists the gain 
in phase across the trench, the gain in Helmholtz free energy, the barrier height (i.e. the gain in Gibbs free energy), and the 
amount of phase that would effectively be generated at the end of the phase slip event (i.e. upon completion of a closed loop, 
the amount of phase generated is the sum of the effective phases) . 



top of the barrier has been passed and the order param- 
eter relaxes to the target metastable state. (To make the 
graph more legible, we have omitted drawing the edge 
that corresponds to this relaxation process.) To find the 
Gibbs free-energy difference between a metastable state 
and a saddle-point state, we need to know the phase dif- 
ference across the center of the trench. To resolve the 
ambiguity of 27r in the definition of G, the phase differ- 
ence can be found by following the wire with no phase 
slip. To further improve the legibility of Fig. ^\ the free- 
energy barriers are listed in a separate table to the right. 
Note, that a phase slip on just one of the wires, being only 
half of the complete process, can be regarded to a gain 
in phase of ±7r for the purposes of calculating voltage, as 
indicated in both the graph and the table. 

Once the table of barrier heights has been computed, 
we can construct a Markov chain on a directed graph, 
where the metastable states are the vertices — in effect, 
an explicit version of our diagram. In general, each pair 
of neighboring metastable states, s n and s n _|_i, are con- 
nected by four directed edges: 



Lop 



"> S71+I 



bottom 



top 



S n +1 S n <- 



bottom 



■+ s n+1 (42a) 
s n +i (42b) 



where the probability to pass along a particular edge is 
given by P(-) = exp— 0AG(a, in which AG(.) may be 
read off from the table in Fig. [§] 

We denote the occupation probability of the n th 
metastable state by o„, where n corresponds to the n 
in the phase constraint @. o n may be computed in the 
standard way, by diagonalizing the matrix representing 
the Markov chain Jig. Each move in the Markov chain 
can be associated with a gain in phase across the device 
of ±7r, as specified in Fig. EI Thus, we may compute the 
rate of phase-gain, and hence the voltage: 



V = — -^—(P(s n ► s m ) - P(s n ► s m j), 

4e f-^ g n , m v top bot " 

{nm) 

(43) 

where the rate prefactor is to be determined, (nm) 
indicates that the sum runs over neighboring states only, 
and g n .m keeps track of the sign of the phase-gain for 
reverse phase-slips: 



9* 



1, if to > n, 
T, if to <C 71. 



(44) 



For the case 6 = 0, and keeping the bottom three states 
only, the voltage generated via sequential phase slips 



turns out to be 

Va=o,seq = — Oe ^\ 3 + -) sinh(//2/ ). (45) 

Having dealt with the case of 8 = (and hence ob- 
tained the value of the resistance at magnetic fields cor- 
responding to resistance minima), we now turn to the 
case of 8 = tt/2, i.e., resistance maxima. 

In this half-flux quantum situation, there are two de- 
generate lowest-energy states, with opposite circulating 
currents. These states are connected by saddle-point 
states in which a phase-slip is occurring on either the top 
or bottom wire. The diagram of the degenerate ground 
states and the saddle-point states connecting them is 
shown in Fig. 1101 By comparing the diagram with the as- 
sociated Table, it is easy to see that the free-energy barri- 
ers are biased by the current, making clockwise traversals 
of Fig. more probable than counter-clockwise traver- 
sals. As there are only two metastable states being con- 
sidered, and as they are degenerate, it is unnecessary to 
go through the Markov chain calculation; clearly, the two 
states each have a population of 1/2. The voltage being 
generated by the sequential phase-slip is then given by 

V s=K/2 , soq = y e Ve- 0e (^-^) sinh(//27 ). (46) 

Vs=it/2, scq is larger than the sum of Vs=o, S cq and 
V£=o, par, so, as expected, the resistance is highest at 
magnetic fields corresponding to 8 — tt/2. For very 
long wires, the perturbation of one wire when a phase 
slip occurs in the other is very small, and therefore we 
expect that the dependence of resistance on magnetic 
field will decrease with wire length. Indeed, for very long 
wires, the difference in barrier heights to sequential phase 
slips between the (5 = and 8 = tt/2 cases disappears 
(i.e. Eq. (HSJ and Eq. agree when b > 1). 



3. Numerical treatment for intermediate-length nanowires 

Instead of using the long-wire approximation, 
Eqs. I|34al34d|) . we can use the exact functions for the 
end-to-end phase accumulation along a wire 6(J(k)), 
and the Helmholtz free energy -F m / sp ( J(k)). By drop- 
ping the long-wire approximation, as the temperature 
approaches T c and the coherence length decreases the 
picture correctly passes to the Josephson limit. In this 
approach, the total current and the phase constraints 
must be solved numerically, as 0(J(k)) is a relatively 
complicated function. Figure IT21 provides an illustration 
of how, for a single wire, the function J{9) depends 
on its length. We shall, however, continue to use the 
barrier-crossing approximation. Because the barriers 
get shallower near T c , our results will become unreliable 
(and, indeed, incorrect) there. 

The form of the order parameter that satisfies the 
Ginzburg-Landau equation inside the wire is expressed 
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in Eqs. l|C20bl IC6(1 . Therefore, to construct the func- 
tions 6{J) and F m/sp (J) [i.e. Eqs. ifCTTl IU2ljl ] . we need 
to find uq(J), i.e., the squared amplitude of the order 
parameter in the middle of the wire. Hence, we need to 
ascertain suitable boundary conditions obeyed by the or- 
der parameter at the ends of the wire. For thin wires, a 
reasonable hypothesis is that the amplitude of the order 
parameter at the ends of the wire matches the amplitude 
in the leads: 

f(z - ±b/2) 2 - g c film( T )^film( T ) ( a 7] 
c wire V / >wirc V / 

For wires made out of superconducting material the same 
as (or weaker than) the leads, this ratio is always larger 
than unity [23| . 

Once we have computed the functions 0(J) and 
Fm/sp{J) for both saddle-point and metastable states on 
a single wire, we can use the phase and total current con- 
straints to build the saddle-point and metastable states 
for the two-wire device. We proceed as before, by con- 
structing a Markov chain for the state of the device, ex- 
cept that now we include in the graph all metastable 
states of the device. By diagonalizing the Markov chain, 
we find the populations of the various metastable states 
and, hence, the rate of gain of O. 

We plot the typical magnetoresistance curves for vari- 
ous temperatures, obtained numerically, as well as dV/dl 
vs. T for various magnetic fields and total currents. No- 
tice that the resistance at 8 = and large total currents 
can exceed that at 8 = tt/2 with low total-current. 



V. CONNECTIONS WITH EXPERIMENT 

In this section a connection is made between our 
calculations and our experiments First, the pre- 

dicted period of the magnetoresistance oscillations is 
compared to the experimentally obtained one. Then, 
the experimentally-obtained resistance vs. temperature 
curves are fitted using our extension of the IZAH Joseph- 
son junction model (for shorter wires) and our extension 
of the LAMH wire model (for longer wires). 



A. Device fabrication 

Four different devices were successfully fabricated and 
measured. The devices were fabricated by suspending 
DNA molecules across a trench and then sputter coating 
them with the superconducting alloy of MoGe. The leads 
were formed in the same sputter-coating step, ensuring 
seamless contact between leads and the wires. Next, 
the leads were truncated lithographically to the desired 
width. In the case of device 930-1, after being measured 
once, its leads were further narrowed using focused ion 
beam milling, and the device was remeasured. For fur- 
ther details of the experimental procedure see Ref. [f| . 
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FIG. 10: Diagram and corresponding Table for the case of magnetoresistance maxima (i.e. 8 — it/2, 3tt/2, . . .). See the caption 
to Fig. |5] for explanation of the diagram. 




«(0) 



FIG. 11: Squared amplitude u(b/2) of the order parameter at 
the end of a wire, as a function of its value uo = w(0) at the 
mid-point of the of the wire, computed using the JacobiSN 
function [see Eq. !C20bH . for the case b = 16. The black line 
corresponds to trajectories that do not go through a pole; the 
gray line corresponds to trajectories that do pass through at 
least one pole. The intersection of the dashed and black lines 
represents those trajectories that satisfy the boundary con- 
dition u(±6/2) = 1. (The intersection of the dashed-dotted 
and black lines represent trajectories that start and stop at 
the same point, i.e., u(b/2) = Uo). 



B. Comparison between theory and experiment 

1. Oscillation period 

The magnetoresistance periods obtained for four differ- 
ent samples are summarized in Table [I] The correspond- 
ing theoretical periods were calculated using Eq. 
based on the geometry of the samples which was obtained 
via scanning electron microscopy. To test the theoreti- 
cal model, the leads of one sample, sample 930-1, were 
narrowed using a focused ion beam mill, and the mag- 
netoresistance of the sample was remeasured. The the- 
oretically predicted periods all coincide quite well with 
the measured values, except for sample 219-4, which was 
found to have a "+" shaped notch in one of the leads 
(which was not accounted for in calculating the period). 
The notch effectively makes that lead significantly nar- 




FIG. 12: Current (in units of the critical current) vs. end-to- 
end phase accumulation for superconducting wires of various 
lengths: 0£ (solid line), 1.88£, 5.96£, 14.4£ (dotted line). The 
transition from LAMH to Josephson junction behavior is ev- 
ident from the loss of multivaluedness of the current, as the 
wire length is reduced. 



rower, thus increasing the magnetoresistance period, and 
this qualitatively accounts for the discrepancy. 

For all samples, when the leads are driven into the vor- 
tex state, the magnetoresistance period becomes much 
longer, approaching the Aharonov-Bohm value for high 
fields. This is consistent with the theoretical prediction 
that the period is then given by Eq. ©, but with I re- 
placed by the field-dependent inter-vortex spacing r. 



2. Oscillation amplitude 

We have made qualitative and quantitative estimates 
of the resistance of two-bridge devices in several limit- 
ing cases. For devices containing extremely short wires 
[b « £(T)], such as sample 219-4, the superconduct- 
ing wires cannot support multiple metastable states, 
and thus they operate essentially in the Josephson junc- 
tion limit, but with the junction critical current being 
a function of temperature given by LAMH theory as 
4(T) = J c (0) (1 - T/T c f/ 2 . A summary of fits to the 
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Sample 


b (nm) 


2a (nm) 


21 (nm) 


Theoretical period (uT) Measured period (//T) 


error 


205-4 


121 


266 


11267 


929.21 


947. 


1.9% 


219-4 


137 


594 


12062 


388.73 


456.6 


12.8% 


930-1 


141 


2453 


14480 


78.41 


77.5 


-1.2% 


930-1 (shaved) 


141 


2453 


8930 


127.14 


128.3 


0.9% 


205-2 


134 


4046 


14521 


47.41 


48.9 


3.0% 



TABLE I: Comparison between measured and theoretical magnetoresistance periods. The geometries of the samples were 
obtained via scanning electron microscopy and used to compute the periods theoretically; see the text for additional details. 



data for this sample, using the Josephson junction limit, 
is shown in Fig. [7] On the other hand, for longer wires it 
is essential to take into account the multiple metastable 
states, as is the case for sample 930-1, which has wires of 
intermediate length. A summary of numerical fits for this 
sample is shown in Fig. 1131 In all cases, only the two low 
total-current magnetoresistance curves were fitted. By 
using the extracted fit parameters, the high total-current 
magnetoresistance curves were calculated, with their fit 
to the data serving as a self-consistency check. As can be 
seen from the fits, our model is consistent with the data 
over a wide range of temperatures and resistances. We 
remark, however, that the coherence length required to 
fit the data is somewhat larger than expected for MoGe. 

VI. CONCLUDING REMARKS 

The behavior of mesoscale NQUIDs composed of two 
superconducing leads connected by a pair of supercon- 
ducting nanowires has been investigated. Magnetoresis- 
tance measurements [f| have revealed strong oscillations 
in the resistance as a function of magnetic field, and these 
were found to have anomalously short periods. The pe- 
riod has been shown to originate in the gradients in the 
phase of the superconducting order parameter associated 
with screening currents generated by the applied mag- 
netic field. The periods for five distinct devices were 
calculated, based on their geometry, and were found to 
fit very well with the experimental results 

The amplitude of the magnetoresistance has been esti- 
mated via extensions, to the setting of parallel supercon- 



ducting wires, of the IZAH theory of intrinsic resistive 
fluctuations in a current-biased Josephson junction for 
the case of short wires and the LAMH theory of intrinsic 
resistive fluctuations in superconducting wires for pairs 
of long wires. In both cases, to make the extensions, it 
was necessary to take into account the inter-wire cou- 
pling mediated through the leads. For sufficiently long 
wires, it was found that multiple metastable states, corre- 
sponding to different winding numbers of the phase of the 
order parameter around the AB contour, can exist and 
need to be considered. Accurate fits have been made to 
the resistance vs. temperature data at various magnetic 
fields and for several devices by suitably tuning the crit- 
ical temperatures, zero-temperature coherence lengths, 
and normal-state resistances of the nanowires. 

As these device are sensitive to the spatial variations in 
the phase of the order parameter in the leads, they may 
have applications as superconducting phase gradiome- 
ters. Such applications may include the sensing of the 
presence in the leads of vortices or of supercurrents flow- 
ing perpendicular to lead edges. 
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APPENDIX A: MULTI-WIRE DEVICES 



In this appendix we give an example of how to extend the theory Presented in this Paper to the case of devices 
comprising more than two wires. In our example, we consider an array of n identical short wires (i.e. wires in the 
Josephson junction limit) spaced at regular intervals. We continue to work at a fixed total current and to ignore 
charging effects. The end-to-end phase accumulations along the wires are related to each other as 

e 2 = e 1 + 26, 
e»3 = e 1 + AS, 

(Al) 



9 n = 9 1 + 2(n - 1)6, 
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T (K) T (K) 

(LHS) Ri = 2882.9 fi, R 2 = 2941.70, f i = 17.3 nA, f 2 = 8.7nA, T cl = 3.147K, T c2 = 3.716K. 
(RHS) Ri = 2912 R 2 = 2912 O, £ i = 10 nA, £ 02 = 9 nA, T cl = 3K, T c2 = 3.65 K. 



FIG. 13: Sample 930-1: Resistance vs. temperature curves. Experimental data (full lines) and theoretical fits using the LAMH 
type model in the intermediate regime (dashed lines). Theoretical curves terminate when the short- wire regime is reached, i.e., 
5£(T) ~ b. (1) Zero magnetic field and low total current. (2) Magnetic field set to maximize the magnetoresistance and low 
total current. (3) Zero magnetic field and 80 nA total current. (4) Magnetic field set to maximize the resistance and 80 nA 
total current. The fit on the left was optimized numerically, and the one on the right was obtained by hand, showing that a 
more realistic value of £oi is remains reasonably consistent. 




FIG. 14: Effective single junction critical current for a multi-junction array, as a function of S. The critical current has been 
rescaled so that J C (S = 0) = 1. Note the similarity with a multi-slit interference pattern. 



i.e., 6 n — 0i — 2(n — 1)6 (for n = 2, . . . , N), where S is the phase accumulation in one of the leads between each pair 
of adjacent wires. The Gibbs free energy of the multi-wire subsystem is given by 

G{I, Qi) = - A J2 cos (0i + 2(m - 1)5) + IO^J , (A2) 

where I is the total current and we are assuming that the wires have identical critical currents. As for the two-junction 
case, this junction array, is equivalent to a single effective junction. Figure^Jshows the critical current of this effective 
junction as a function of d for devices comprising 2, 5, and 15 wires. The magnetoresistance of such a device then 
follows from IZAH theory, i.e., Eq. lpTT|) . 



APPENDIX B: PHYSICAL SCALES 



It is convenient to express the results of the long- wire model, Eqs. I|41I45I46|I . in terms of macroscopic physical 
parameters. Following Tinkham and Lau Il9l | , we express the condensation energy scale per coherence length of wire 
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as 

'■"■w-^iliRr < B1) 

where t = T/T c , i?N is the normal-state resistance of the device, and R q = h/4e 2 ss 6.5 kf2 is the quantum of resistance. 
The LAMH prefactor for sequential phase slips then becomes 

bVT—t ( 8V2£ \ 1/2 8k B (T c -T) 
n -W^)[Wc) ^ ' ( } 

and for parallel phase slips becomes 

The remaining parameters in the model are Rn, T c and £(T = 0). The normal-state resistance and the critical 
temperature may be obtained from the R vs. T curve. The coherence length may be obtained by comparing £ 
obtained from the critical current at low temperature, via 

- 2 W7rE (B4) 



c 3V3 $o ' 

with £ obtained via Eq. (|B1(1 . 

In experiment, it is expected that the two wires are not identical. The long-wire model can be easily extended 
to this case. The number of parameters to be fitted would then expand to include the normal-state resistance for 
each wire (only one of which is free, as the pair are constrained by the normal-state resistance of the entire device, 
which can be extracted from the R vs. T curve), a zero-temperature coherence length for each wire, and a critical 
temperatures for each wire. 



APPENDIX C: LAMH THEORY FOR A SINGLE BRIDGE 

In this appendix we reproduce useful formulas from LA and rewrite them in a way that is convenient for 
further calculations, especially for numerical implementation. As in the case of single-wire LAMH theory, one starts 
with the Ginzburg-Landau free energy 

F = f 2 c#| 2 + f + iw| 2 dz. (CI) 

J -b/2 2 2m 

The relationships between the parameters of the Ginzburg-Landau free energy (a and /?), coherence length £, the 
condensation energy per unit coherence length £, the critical field H c and the cross-sectional area of the wire a are 
given by ^- = = £/£ and £ 2 = 2 m\a I ' F°ll° w i n g McCumber ^tJi it is convenient to work in terms of the 

dimensionless units obtained using the transformations: \ip\ 2 — > ^|V ; | 2 7 z y 2 "1" z, and b ^ b/{ = y^Pfc. In 
terms of these units, the free energy becomes 

F = 2£ Q(l - IVI 2 ) 2 + |VV-| 2 ) dz. (C2) 

The Ginzburg-Landau equation is obtained by varying the free energy: 

SF = -^0 + |V'| 2 ^- V 2 ?/> = 0. (C3) 

By writing ip = fe 1 ^ and taking the real and imaginary parts of the Ginzburg-Landau equation one obtains 

-/ + f + (4>'?f = /", (C4) 
20'/' + 4>"f = 0. (C5) 
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FIG. 15: "Mechanical potential" U[u — f 2 ] at an intermediate value of the dimensionless current, plotted as a function of 
amplitude squared to make comparison with Fig. Illl more convenient. 



From Eq. I|C5|) . one finds the current conservation law: 

/V = J, (06) 

where J is identified with the dimensionless current -^(^Vtp — tpWifj*)- The physical current (in stat-amps) is given 
by / — JcH 2 a£_/^>o- Expressing <p' in terms of J, Eq. I|C4|) becomes 

f" = -f + f + ^ = ~U(f), (C7) 
where the effective potential U(f) is given by 

72 f2 fi 

u( f) = w + L-L. (08) 

Following LA, Eq. I|C8|I can usefully be regarded as the equation of motion for a particle with position f(z), where z 
plays the role of time, moving in the potential U(f) [l4|. Before proceeding to find the solution of this equation, we 
pause to consider the type of trajectories that are possible. Later, it will be demonstrated that at the edge of the wire 
/(±6) > 1, so the particle starts to the right of the hump; see Fig. El If the total energy of the particle is less than 
the height of the hump, the particle will be reflected by the hump. If, however, the particle starts with more energy 
than the height of the hump, it will pass over the hump and be reflected by the J 2 /2f 2 dominated part of U(f). 

The equation of motion can be solved via the first integral (i.e. multiplying both sides by / and integrating with 
respect to /): 

(/' 



E="J- + U(f) => f' = y/2(E-U(f)), (C9) 
where E is a constant of integration (i.e. the energy of the particle in the mechanical analogy), which gives 

df f f W (cio) 



/ 

/o ~^J2(E -[/(/)) J fo ^2/ 2 £- J 2 - J 4 + / 6 /2" 

It is convenient to apply "initial" conditions at the middle of the wire, where f(z = 0)=fo, and to integrate towards 
the edges. We require that the particle come back to its starting point after a "time" b, i.e. at the edges of the wire 
the amplitude of the order parameter must match the boundary condition. Therefore, the middle of the wire must 
be the turning-point for the particle, i.e., at z = we have E — U(fo). 

What follows next is a series of manipulations via which one can express solution for f(z) in terms of special 
functions. 

Step 1: substitution: f 2 ^u 

z =ir - (en) 

2 J U0 y/2Eu -J 2 -u 2 + u 3 /2 
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FIG. 16: Squared amplitude u of the order parameter as a function of position along the wire for the two types of solution: 
metastable and saddle point. 



Step 2: substitution: u — > uq + e 




(C12) 



a 



-'•-=/ 37^7^/2 ( C13 ) 



(e(e + /? e + e72))V2 



(e(e + 13 + v//32-2a)(e + /? - v 7 /? 2 - 2a)) 1 / 2 



-Ml 

%/2de 



(e(e-u 1 )(e-u 2 ))V2 



(C14) 



(C15) 



Step 3: substitution: e — > u\z 2 



_ 2V2 rv -=r- 



2 V / 2 T , 11 . . 



EUipticFrArcSinr A /^-^l,— 1 (C17) 

/U2 V Ul U2 

The following definitions have been used: 

3 

a[u ] = J 2 /u - u + u , P[u ] = -uo - 1, (C18) 



Ui[a,/3] = -(3 - v//3 2 - 2a, «2 [">/?] = -/3 + V/? 2 - 2a. (C19) 

By inverting relation (|C17J| one obtains an explicit equation for the amplitude of the order parameter as a function 
of position along the wire (see Fig. 116(1 : 

f 2 (z) = u a + u\ sin 2 TjacobiAmplitudefzw — 2 -, —11 (C20a) 

V 2 U2 

= u + iiiJacobiSn 2 [zJ — .— 1 (C20b) 

V 2 U2 

The end-to-end phase difference along the wire may be found by using the current conservation law. Thus one 
obtains 

6/2 J f b/2 dz 

dz = 2J — n 2r =^-t - (C21) 



b/2 P{z) Jo u + uiJacobiSn 2 [ Zy /%, £ 
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The Helmholtz free energy can be found by substituting the expressions for f(z) and 4>'(z) into the expression for 
the free energy. One then obtains 

,6/2 i 

F = 4£ dz(- -2f + f 4 + J 2 /u + u -u 2 /2), (C22) 
Jo * 

where E was expressed in terms of uq. Eqs. (JC21I IU22|l provide expressions for 8 and F which are true regardless 
of the length of the wire, and therefore may be used as a starting point for computing the Gibbs free-energy of the 
various metastable states subject to the total-current and t he phase constraints. 
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obvious whether ±(6/2)-^/it2/2 lies in the first period, as 
can be seen from Fig. 1111 As the trajectory must be sim- 
ply periodic, zyf U2/2 must intersect either a zero or a 
pole in the first unit quarter cell of the JacobiSn func- 
tion. Now, we are only interested in trajectories that es- 
cape to / — > 00 [as /(±6/2) is assumed to be greater 
than or equal to unity], so a pole must be intersected. 
(However, being outside the first period is unphysical, as 
it means that somewhere along the wire / = 00.) There 
are exactly two poles in the first unit quarter cell. They 
are located at 2vi +V2 and V2, where vi = ~K{ui/ui) and 
V2 = iK(l — W1/U2), in which K(-) is the complete elliptic 
integral. So, instead of checking whether ±(6/2)^/112/2 is 
outside the unit quarter cell, we can just determine which 
pole z\J U2I2 intersects and then see if ±(6/2)-^/ U2/2 lies 
beyond that pole or not. 



